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We develop a 3+1+1 covariant formalism with cosmological and astrophysical applications. First 
we give the evolution and constraint equations both on the brane and off-brane in terms of 3- 
space covariant kinematical, gravito-electro-magnetic (Weyl) and matter variables. We discuss the 
junction conditions across the brane in terms of the new variables. Then we establish a closure 
condition for the equations on the brane. We also establish the connection of this formalism with 
isotropic and anisotropic cosmological brane-worlds. Finally we derive a new brane solution in the 

. . . framework of our formalism: the tidal charged Taub-NUT-(A)dS brane, which obeys the closure 

^— ' condition. 

^ I 

O '• PACS numbers: 

(N ■ 

. I. INTRODUCTION 



In recent years the idea of exploring the possibility of the gravitational interaction acting in more than 4-diniensions 
[ (4d) attracted a lot of attention. In particular, one of the simplest such models arises from the curved generalization 
of the one-brane Randall-Sundrum model Q (for a review see [^), where gravity acts in 5-dimensions (5d), however 
p— i' standard model fields are confined to the 4d brane. Nonstandard model fields can occur in the 5d space-time. A 
O 2^2-symmetric embedding looks natural only if the brane is envisaged as a boundary' otherwise generic asymmetric 
^ embeddings should be allowed. The generic dynamics was given in Refs. Q and Q in a 5d-covariant approach. 
5-H Although promising at the level of allowing new degrees of freedom, which seem to be adjustable to explain for 
I example dark matter the complexity of the brane- world dynamics also represents a major impediment in 

] obtaining simple exact solutions or in monitoring the evolution of perturbations. Therefore new approaches to 
CNJ ■ gravitational dynamics in brane-worlds are worth to study. 

K>" ] In Ref. 0] the 5d space-time was foliated first by a family of 4d space-times and a second family of 4d space-like 
i: hypersurfaces (see Fig 1. in Ref. @). Such an 5d space-time is double foUable. This structure of the 5d space-time 
allowed to describe the gravitational degrees of freedom in terms of tensorial, vectorial and scalar quantities with 
respect to the 3-space emerging as the intersection of the two hypersurfaces. They represent the gravitons, a gravi- 
photon, and a gravi-scalar and are given by quantities with well-defined geometrical meaning, namely the tensorial, 
vectorial and scalar projections of the spatial 4d metric and their canonically conjugated momenta: the extrinsic 
curvature (second fundamental form of the 3-spaces with respect to the temporal normal), normal fundamental form 
' and normal fundamental scalar. The evolution equations for these variables were given a la ADM both on the brane 
and outside it. An extension of this formalism towards a full Hamiltonian description was advanced in Ref. [lo| . It 
has been shown, that among all gravitational variables on the brane, only the momentum of the gravi-vector has a 
jump across the brane, related to the energy transport (heat flow) on the brane. 
}^ • The formalism presented in Refs. however is not straightforward to be applied in brane cosmology. A formal 

5^ I difference would be the definition of the time derivative. In Refs. this was defined in the tradition of canonical 

gravity as a Lie-derivative taken along the temporal direction (which is not necessarily orthogonal to the 3-space) 
projected to the 3-space, whereas in cosmology by tradition the time derivative is defined as a covariant derivative 
along the normal to the 3-surfaces (this derivative happens to enter only in expressions projected onto the 3-space). 
Obviously this mismatch in the definitions of the time derivatives is not crucial, as the two definitions differ only in 
terms taken on the 3-space. However the formalism developed in Refs. 0] and relying on the double foliability of 
the 5d space-time is unable to deal with the possible vorticity of the word-lines of observers. 

In this paper we develop a formalism overcoming this inconvenience and derive the full set of evolution and constraint 
equations governing gravitational dynamics in a 3+1+1 covariant form. Provided the space-like normal n has vanishing 
vorticity on a hypersurface (the time evolution vector can still have vorticity), the formalism becomes suitable for 
describing gravitational dynamics on the brane (then n becomes the brane normal). In this sense the formalism is a 
generalization of the brane 3+1 covariant cosmology [ll| (which in turn is a generalization of the general relativistic 
3+1 covariant cosmology 12]-il5'|). 

This newly developed formalism also generalizes the s+1+1 covariant brane- world dynamics developed in Ref. [^. 
Both the vector field n and the time evolution vector u are allowed to have vorticity in the present formalism. Though 
observational evidence suggest that the directly detectable 3+1 part of the universe is best described by a Friedmann 
brane with perfect fluid, when discussing cosmological perturbations, the vorticity of u should be allowed, similarly 
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as in existing formalisnis of covariant cosmology in both general relativity and brane-worlds. We can also argue for 
the need of keeping the vorticity of n. One reason would be, that if it is vanishing at some initial instant, it should 
stay zero; and the formalism should be able to handle its "evolution". Secondly, and more important, the vorticity of 
n should not necessarily vanish at other locations, than the brane, a gauge freedom worth to explore. 

We note that a lower-dimensional, 2+1+1 formalism was developed in Refs. [l^, [l^ and applied in the general 
relativistic covariant perturbations of Schwarzschild black holes and rotationally symmetric space-time; then for 
investigating spherically symmetric static solutions in f{R) gravity theories in Ref. [18|. 

In general relativity the important topic of cosniological perturbations, related to both the Cosmic Microwave 
Background and structure formation, has a rich literature, from which (without claiming completeness) we mention 
Refs. [l^-[23], all based on the 3+1 covariant approach. 

Brane-world cosmological perturbations are equally important, however additional difficulties emerge due to the 
complexity of brane-world theory and the impediment to predict and perform observations on the brane. At a technical 
level, the latter is obstructed by the lack of closure of the perturbation equations on the brane. Despite impressive 
developments [l^l-jiOlj many questions remain unanswered. Although we cannot overcome well-known difficulties, we 
expect our new formalism will provide the most convenient and complete tool-chest for approaching the problems. 
We also foresee the possibility of important applications for brane-world exact solutions. 

We establish the generic 3+1+1 covariant formalism in Sec. |lTl by defining the kinematical quantities and the 
decomposition of the Weyl and energy-momentum tensors. We also relate the curvature and Ricci tensors and the 
3-dimensional (3d) scalar curvature to kinematic, non-local (Weyl) and matter-defined variables. In Appendix VK\ the 
commutation relations among all types of derivatives emerging in the formalism are given. 

Sec. mil contain the full 3+1+1 decomposed covariant gravitational dynamics and constraints, together with the 
available matter field evolutions. In Appendix |B] we discuss the gauge freedom in the frame choice and give the 
transformations of all relevant quantities under infinitesimal frame transformations. 

Taking into account that the brane is a 4d time-like hypersurface, in Sec. IIVI we discuss the decomposition of the 
Lanczos equation and of the sources of the effective Einstein equation. Then we specify the generic evolution and 
constraint equations on the brane, expressed in terms of quantities defined on the brane. These equations arise from 
combinations of the equations given in Section [III[ evaluated on the brane. Appendix [C] contains the brane equations 
for an asymmetric embedding. We continue Section HVl bv specializing to a symmetrically embedded brane, by taking 
into account the Lanczos equation. Then we conclude the section by deriving a closure condition for the equations 
on the brane. 

Sec. |V] contains the derivation of the most important cosmological equations to a hypersurface (a generic brane), 
then the discussion of the particular case of cosmological symmetries and perfect fluid source on the brane. As a test 
of our formalism we recover the Friedmann, Raychaudhuri and energy balance equations and compare them with the 
relevant results of Ref. 3- In subsection IV Bl and Appendix [P] we also relate our formalism to anisotropic brane-world 
cosmologies. 

Sec. IVII contains an astrophysical application of our formalism devoted to locally rotationally symmetric (LRS) 
space-times, at the end of which we recover a new brane solution with NUT charge. This solution obeys the previously 
derived closure condition. 

Sec. IVIII contains the Concluding Remarks. 

Notations. Quantities defined on the 3-space orthogonal to both u° and carry no distinguishing mark and 
the 3d metric is denoted hab- Quantities defined on the brane carry the pre- index the only exception being the 
4d metric gab- Quantities defined on the full 5-space carry a distinguishing ~ mark. Exceptionally, other quantities 
also carry the distinguishing~mark. These are (a) the 3+1+1 decomposed components of the 5d energy-momentum 
tensor, which are defined on the 3-space orthogonal to both m° and n°, and (b) the kinematic and extrinsic curvature 
type quantities related to another singled out spatial vector e° in Sec. IVII Calligraphic symbols denote 3+1+1 
decomposed components of the 5d Weyl tensor. Whenever possible, identical symbols are used for quantities related 
to the temporal vector field u"" and the brane normal n°' , the latter distinguished by a^mark. We denote the average 
of a quantity / taken over the two sides of the brane by (/), and its jump by A/. Angular brackets ( ) on abstract 
indices denote tensors which are projected in all indices with the metric hab, symmetrized and trace-free. Round 
brackets ( ) and square brackets [ ] on indices denote the symmetric and antisymmetric parts, respectively. 

II. THE 3+1+1 COVARIANT FORMALISM 

A. Decomposition of the metric 

Let M° = dx"" /dr and n° = dx'^/dy be a time-like and a space-like vector field in the 5d space-time, with r 
and y the affine parameters of the respective non-null integral curves. They obey the normalization conditions 
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FIG. 1: Elements of the 3+1+1 decomposition of the 5d space-time geometry with metric Ijab and compatible connection V. 
The 4d brane with normal n" has induced metric gat ~ Ijab — naUt, and extrinsic curvature '*'i('ab = glgt^ end- The 4-velocity 
field u° of observers in the brane defines local 3d orthogonal spatial patches, hypersurface-forming only when the vorticity 
^ab ~ 0. The other kinematical characteristics of w" are the expansion and shear aab- The expansion, shear and vorticity of 
n" are G, Gab and tDab, the latter vanishing on the brane. The temporal, off-brane and 3d covariant derivative are shown for 
the vorticity 1-form uja- 



~u°'Ua = 1 = n°'na and the perpendicularity condition u°"na = 0. The 5d metric is decomposed as 

gab = naUb + gab , 

with 



gab 



-UaUb 



(1) 



(2) 



the metric on the 4d temporal leaves y =const (with the brane at y = 0) and the spatial part hab of this metric 
obeying u°-hab — n'^hab — 0. We denote by Sabc the volume element associated with hab- The 4-vector u°- represents 
the time-like velocity field of brane observers (see Fig 1). 

We employ three type of derivatives, all associated with projections of the 5d connection V^. A dot and a prime 
denote covariant derivatives along the integral curves of m° and n", respectively, while D is the 3d covariant derivative 
compatible with the metric h: 



fb..c = WVaTb^c , 

TLc = n^^aTb^c , 

DaTb..c = ha^hb\-K'^M„j 



(3) 
(4) 
(5) 



Note that the Z?-derivative is the same as introduced in general relativity employing the corresponding projection 
of the V-derivative (V being the connection compatible with gab)- This is, because both generate the covariant 
derivatives formed with the connection compatible with hab- Concerning the time-derivative defined above, except for 
scalars, it differs from the corresponding general relativistic time derivative employed in the 3+1 covariant formalism 



(which is defined with V) 
definitions agree. 



in n°- and u°- terms. Nevertheless, when projected to the 3-manifold with the two 
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B. Kinematic quantities 

We introduce the kinematic quantities through the decomposition of the 5d covariant derivative of n° as 

VaUb = -UaAb + KUaUb + KUaUb + UaKb + LaUb + Kab , (6a) 
VaTlb = UaAb + KUaUb + KUaUb ~ UaKb + LaUb + Kab , (6b) 

where 

A-a = U{a) , Aa^ u'^a) ' 
Ka = "(a) , Ka^n^a) , 

K = n''u[ , K = u^ib , 
La = h^'^n'^VcUd , 

^i^ab = DaUb , A'ab = DaUb . (7) 

As a rule, an overhat is used for kinematical quantities related to the vector field n° in order to distinguish them 
from the similar kinematical quantities related to the vector field w". Here Aa is the acceleration. All scalars, vectors 
and tensors in the above decomposition are defined on the 3d manifold orthogonal both to n" and The tensorial 
expressions Kab f^nd Kab can be further decomposed into (trace-, trace-free symmetric and antisymmetric) irreducible 
parts as follows: 

Kab = -^hab + CTab + ^ab , (8) 
Kab = ^hab+^ab+UJab , (9) 

where we have defined the expansion, vorticity and shear of the vector fields and as 

e = W^Ua , UJab = D[aUb] , (Tab = D^a^b) , 

e = W^Ua , UJab = D[anb] , O^ab = D (^a^b) ■ (10) 

The antisymmetric 3d tensors ujab and ujab can be also encoded in the vorticity vectors to"' — e'^^'^ujbc/'^ and cD" ~ 
e°''^Wbc/2. When the vorticities of both vector field and u°- vanish Wab ~ ^ab — 0, the tensorial expressions Kab and 
Kab are symmetric and they represent the two extrinsic curvatures of a 3d hypersurface. The condition Qab = is also 
necessary in order to have a brane at y = 0, but not sufficient. Indeed the brane is a 3-1-1 dimensional hypersurface 
which can exist only if the higher-dimensional vorticity of its normal ujat = ^[ai^b] vanishes. This condition translates 

into = '^<7f,]'^^c'^d = —U[a (^Kb] + Lb]j +ujab, therefore (due to Frobenius' theorem) the necessary and sufficient 
conditions for the existence of the 3-1-1 brane are 



^ab\y=o = 



(11) 



In summary, the independent components of VaUfc are expressed by three scalars (K, K, O), four 3-vectors 
{Aa, Ka, La, uja) and a symmetric trace-free 3-tensor (iJafc). The corresponding decomposition of Va?i6 consists 
of the three scalars {K, K, Q), four 3-vectors {Aa, Ka, La, uja) and a symmetric trace-free 3-tensor {dab)- The 
irreducible decompositions of the covariant derivatives VaWfe and Van^ have therefore 20 independent components 
each (due to the constraints w^'VaWb = and n^'Vanb = 0). 



C. Gravito-electro-magnetic quantities 



The non-local gravitational properties of the 5d space-time are carried by the 5d Weyl tensor, the principal directions 
of which lead to a classification scheme generalizing the general relativistic Petrov classification (41,]. Here we are 
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interested in the 3+1+1 decomposition^ of Cabcd, which can be given in terms of the quantities (with a total of 35 
independent components): 



£ 








= iefe "'Cabcdu'n" , 


J^ki — Cabcdh(^^v!^hi-^'^n'^ , 


£k 


= Cabcdh''n''u-n'' , 


£k ^ Cabcdhk''u'n''u'^ , 


£kl 


= Cabcdh(^^."'n'hi-^''n'^ 


£kl — Cabcdhi^j^'u^hi^'^u'^ , 


Ukl 


abi^ d 
— 2 ') ^abcdn 


, 'Hkl = -£(fc '^^hi-^'^Cabcdu'^ 



(12) 

We note that aU tensorial quantities defined above are trace- free (from the properties of the Weyl tensor), and further 
in the tensors Fki-, Hki and Hki the brackets ( ) are equivalent with the round brackets ( ). The Weyl tensor in terms 
of the quantities defined in (IT^ is 



Cabcd — —2 {£d[ahb]c ~ £c{ahb\d) + 2 ^fd[a^fc]c " £c{ohb\d^ — 2^^c[a^6]d 

+2 {e^d '"H-i\anb\ + ^ab ""Hiicnd^j - 2 (e^^ ''"Hi[aUb] + e^b "Hi[cUd]) 

+ 2 (ncn[a£b]d - ndn[a£b]c^ + 2 {UcU[a£b]d - UdU[a£b]c) 
-2 (Uc"-[a-^fc]d - Ud«[a-^6]c + "-cU[a-^fc]d " "-<iW[a-^6]c) 

- {nycUa\£bdk + U[cnb]£adk) T-L'' - {u[dna]£bck + n[dUb]eack) T-L'' 

-2 {U[cnd]£abk + "[a^bjEcdfe) 'H*' - 2 (£[a^[c"-<i] + '^^[c/irf] [a"b] ) 
-2 {£[ahb][cUd] +£[chd][aUb]j + 4f UfaTlfc] U[cn<j] 

-2 {ncUdn[a£b\ - Ucndn[a£b] + naUbUy^Ed] - Uanbn[c£d]^ 

+ 2 (ucndU[aSb] - ncUdU[a£b] + UaUbU^cEd] - naUbU[c£d\) 

2 2 

--£{ndn[ahb]c~ncn[ahb]d) + ^£{udU[ahb]c-UcU[ahb\d) ■ (13) 

This relation generalizes the general relativistic 3+1 covariant decomposition of the 4d Weyl tensor Cabcd (which 
has only 10 independent components), where only two tensors Eu = Cabcdh^f.°'u''hj^s^'^u'^ and Hki = ^e^^. '^''h^^'^Cabcdu'^ 
appear, the electric and magnetic parts of the Weyl curvature. On the brane the relation between the set of variables 
Ski and Hkl and the variables Eki and Hki is given by 

Sab = Eab + ^Sab ^ ^ + ^ ^ + ^^cia'^b)' + ^^(a^'b) , (14) 

■Hab Hab ~ '"^'^b)cKd ■ (15) 

D. Decomposition of the energy-momentum tensor 

The 5d gravitational dynamics is governed by the Einstein equation 

Tab + TabS {y) , (16) 



Gab = -^9ab + 



where k denotes the 5d coupling constant. The sources of gravity are the 5d cosmological constant A, the regular 
part of the 5d energy-momentum tensor Tab and a distributional energy-momentum tensor with support on the brane: 

Tab = -^9ab + Tab ■ (17) 



^ The case of a generic n + 1 decomposition of the Weyl tensor was discussed in Ref. |42|| . 
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Here A is the constant brane tension and Tab describes standard model matter fields on the brane^, decomposed with 
respect to a brane observer with 4-velocity as 

Tab = PUaUb + q{aUb) + phab + ab ■ (18) 

The quantities p, qa, p and iTab are the energy density, the energy current vector, the isotropic pressure and the 
symmetric trace-free anisotropic pressure tensor of the matter on the brane. 

We decompose the regular part of the 5d energy-momentum tensor relative to and n°- as: 

Tab = PUaUb + Iqf^aUb) + '2qU(^anb) + phab + TTnaUb + 2TT(^anb) + ab , (19) 

where p is the relativistic energy density relative to u", p the isotropic pressure, 5a and g the relativistic momentum 
densities on the 3-space and in the direction n°. The quantities tt, tTq and the symmetric and trace- free tensor iVab are 
related to the scalar, vectorial and tensorial components (with respect to the 3d space) of the 5d anisotropic pressure 
tensor, which is 

3{tt—p) _ p—TT ^ 

riam + 2TT(anb) H ~'^'^b + T^ab ■ (20) 

By employing the definitions given in this Section we give the commutation relations among the temporal, off-brane 
and brane covariant derivatives in Appendix A. 



E. The Gauss equation and its contractions 



We define the local 3d curvature tensor TZabcd of the space orthogonal to both and n° as 



1 



resulting in 



D[aDb]Vc- UJabVi^c) +^abVl^~^ = -TZabcdV^ , 



Tlabcd = h^'h^^hJ'h^^Ryki - (DaUc) {DbUa) + [DaUd) {DbU^) 

+ {DaUc) (DbUd) - {DaUd) (DbUc) ■ 



(21) 



(22) 



This is a natural generalization of the definition used in general relativity [2^ , 14J] , |45| . By the definitions of the 
kinematical, gravito-electro-magnetic and matter variables, we have 



TZa 



bed 



hrtn.h, 



afl^d 



-- (e^-e^j -£ + A + K^{p+p-n) 

-2 {£d[ahb]c - £c[ahb]d) + 2 {£d[ahb]c - £c[ahb]dj ^ {hd[a'n:b]c " hc[a1^b]d) 

29 



[{<^a[c + ^a[() hd]b + {crb[d + hc]a] ~ 2 (a a[c + i^a[c) {o' d]b - ^d]b) 



26 



H ^ [{^a[c + t^a[c) hd]b + (CT6[d + hc]a] + 2 (da[c + '^afc) (f?d]6 - '^d]b) 



(23) 



For DGP / induced gravity type models |43|, Tab should be replaced by Tab — {l/'^'^) Gab, with Gab the Einstein tensor constructed 
from the metric Qab and 7 the dimensionless induced gravity parameter. Randall-Sundrum type brane-worlds are recovered for 7—^0, 
on the e = — 1 branch. 
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and the local 3d Ricci tensor and Ricci scalar become 







~2 


26 




"3" 



bed 

[ d \ / 

+£ac — Sac + "^""ac ^ (o'ac + ^ac) + i'^ ac + ^ac) 

+0'afcO-c^ ^ l^b^^hac + ^a^c - 2o-[^^Wc]{, 

-f?ahf?c ^ + i^bi^^Kc - ^aP^c + 2CT[„ ''t2ic]6 , (24) 

7^ = /^''"7^ac = -26: + A + (p + p - i) - ^ (6^ - §2 

-2tJaW° + CTabO-"^ + 2a}aW° - CTabCT^'' . (25) 

Eq. (|25p can be referred as a generalized Friedmann equation in the 5d space-time, as will become evident in Section 
IVl The general relativistic counterpart is presented in Refs. [ist - 

III. 3+1+1 COVARIANT DYNAMICS 

The full set of the 3+1+1 covariant evolution and constraint equations for the kinematic, gravito-electro-magnetic 
and matter variables are given by the projections of the Bianchi and Ricci identities for and n". These equations 
hold as the main result of the paper, and they are presented in the following three subsections. In particular, the first 
subsection contains all Ricci identities; the second subsection contains twice contracted Bianchi identities, which by 
virtue of the Einstein equations describe evolutions for the energy density and currents; while the third subsection 
contains the rest of independent Bianchi identities. A related Appendix |B] gives the transformation rules under a 
frame change for the totality of the kinematic, gravito-electro-magnetic and matter variables, to linear order accuracy. 

A. Ricci identities 

The Ricci identity for u"^ gives the following independent equations: 

= k + k' + AaA" - k^ + LaK" - kaK"" - Lak" 

+f - I + y (p - i + 3p) , (26) 
02 

= 6 - i:>"Aa + — + 6^ - A-'Aa - 2WaW" + (T^bCT"'' - f 
o 

k-'Ka + K'^K + L^Ka) - | + y(/5 + i+ p) , (27) 

= K\a) - ^(a> +{K+ I I Ka+(K~^]ka+ £a 



3 / V 3 



-2 



-K{Aa+Aa]+ (i^ba + Tba) [K' - K' j - yi, , (28) 

= L(,) + Dak + ( A' + 1) La + (k +^^Aa+(K- 1^ ka 

-2 

+ (Wab + ^ab) A" - {uJab + (Jab) [K" - j + - y , (29) 
- LJ^a) - ""^DcAd + kQa + ^t^a " (TafcCj'' + ^'^ (^K^kd + KcLd + k^Ld) , (30) 



26 

= &(ab) - Di^a^h) + -^(^ab + Kaab- Ai^a^h) + Kl^aLb) 

^ ^ ~2 

+K{aKb) + L{aKb) + W(aWf,) + CF c{a<^b) + ^ab " ~^^a.b ' (31) 

= 6'- D'^Ka -(^K-^% + {K'' + i») A„ - A» - L„) - 2a;„w» + - , (32) 

= L\^^ - D,K+ I) - I^K- i„ + l) ^« 

~2 

+ (?a6 + ^ab) (K'' + L'') - (Wafc + - f „ + y , (33) 



- - O'ab + UJ{a^b) + ^cia^b) " + {/^b)d + (^cia^b)" + ^ab , (36) 
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-(i^'^/f.-ifu,-™,) + | + y + , (45) 



(34) 



= - -^D^Kb - (a, + 1,) {K^ - L,) - - |) 

= <^\ab) - D(aKb) - A<6 {Ka) " ^a)) + ^(6 (-f^a) + ^a)) + f <7a6 



= £fc + 2 + I j Wfc + 2 (^i^ - 1^ Wfe + £fe "^WaWb - Gka^" + <Tfe„w" - 6^ "^06^/ + "^fe , (37) 

= i^"Wa - AaW" + (ifa - ia) , (38) 
= -D(cWfc> + SabikD^CTc)' + 2-4(cWfc) - 2Ki^j::f,) - ^(cWfe) + £(j.°''CTc>fe-f'a + ^fec , (39) 

^ t. 2 26 t 

= D^aab + e„ ""D^ujk - -^D^Q + — L„ - e„ ^^'^L^Wfc 

+2£„ {A^UJk - KcLJk) - Oabiy + <?a + — 96 • (40) 

The Ricci identity for n°' gives the foUowing independent equations: 

= § - D'^k^ + + I j e _ _ ^„ + > (k^ + _ 2a;„a;" + a^bo"^ - K^q , (41) 

= ^ik) - \ek "'Dakb + (lb + Ab) + i„) + ( 



+ + ^ j (^ab + ^(oWfc) - Wc{a^6)'' " <^(a'^'^b>d + '^ciaO^b)" + ^ab , (43) 

~2 

+ (W6a + afaa)(if^-/^^) +4-y§a , (44) 
= e'-£'M„ + ^-e/<: + lM„-2a)„u;" + a„ba'^''+£: 



= - l^a """DcAd - KiOa + ^S, - dab^' - ^s, (^K,Kd + K,Ld + K,L^ , (46) 



~2 
K 



-Kl^aKb) + L{aKb) + i^{a^b) + (fc{aPb) + ^ab + "3"'^a6 ' (47) 
= D'^Ua + - UJ'' (ka + ia) , (48) 



; I 2-29 , 

= D''aab + £a D,u:u~-DaQ + —La-e,'''L,Uk 

-2e„ ^'^ (leWfe - i^cWfc) - CTafci' - fa - • (50) 



B. Conservations laws 

— a — 

The twice-contracted 5d Bianchi identities imply V Tab = 0, which can be decomposed into the projections taken 
with u, n and h, respectively: 

/5 + g' + i:i'^9„ + p(x + e) + x5? + ep + i="''fT„6 

'q{2k {2Aa - Aa) + {La + Ka) , (51) 

= q + lr' + D^na + n (q - k'j - Kp- %p - ^''^dab 

+q{2K + 9) - i'' (2Aa -Aa)-r (ka - La) , (52) 



= ^(fe) + ^'(/c) + ^kP + D°-^ak + ^T^fc + ^'Jfc - k-Kk + Kqk + pAk + nAk 

-p (Ak - Ak) + q-^LOak + g Vafc + ^"-Uak + ^^CTafe +q{kk+ Kk) - ^ ak ( A" - A" ) . (53) 



The first of these equations is the continuity equation, as can be easily verified in the homogeneous, isotropic case 
(q = g° = tt" = 7?°^ = and 9 = iH) and for K = Q. The ensemble of the equations represent an incomplete set of 
evolution equations for the 5d matter (there are no evolution equations for p, tt, tTq, TTab)- 

C. Bianchi identities 

The equations independent from Eqs. ([5T|) - (l55|) arising from the 5d Bianchi identities are: 

= £ - D'^Ea + + - -^afc?"' + iUa^" - £a (2A"^ + L°) - 2£aA'' 

2'k 4k^^ ,„ 21^^ _\ 27? 



3 



eg ~ ^gaA« - {2k'^ + L'') - ^i.^a'^" , (54) 
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4£ ^ 1 ^ ^ ^ ^ ^ 



-El 



afc 



2k 



2 . 



2k 



2k 



, 2kV- e\ _ 2k\^ ^ 
(p + tt) Lfe 3~ \ 3" 1 * 3~ ^ 

2k^^ 2k% 2K^e^ 2K^g^ 

^TTfcaii + (Wfca + CTfeaj TT H ^ — TTfc H ^^4^, 



(55) 



Aa£b — La£b 



— £ 



ab 

k '-'ac^b 
2 



^Sfe + — Q'Wfc + yEfe ""TTacfTfe + — (tT - p) Wfc + y TTfcaO; , 

f (fe) + f - DkS - SkaA°- - EkaA°- + KEk -K£k + ^ (oSk + QSt) 

+2 {oJka + Tka) + 2 + aka) + Tka (i^" + + ^£fea6^" (i^' - i^") 

- ^5 (Afe - Ife) - y Dfc (p - ?f + p) + y -D"^afc 
2k^ ^ \ K^ „ 



(56) 



(57) 



= £(k3) - ^'i 



'{ki) ~ £ab{k 



jja^b 



j) 



-Dtk£. 



K 



e 



£kj + {K + Q)£kj 



+ '^^kj +2\k-'^ Tki - £(k {a^) - \h) + \^{k (4^.) - H - 3Kk^) 

3 3 

+ 7{^{k'^3) + '2^a{k'^ j) + " ^k)aHb + ^ (fe " (wj)a - 3aj)a) - J"^^" (wfc)^ - dk)a) 



2k^ 



At^ - e^^ (/^Tfc + Lb) + ^TT^kj) + ^D^kQj) + ^©^r 



2k^ 



2k^ 



2k^ 



jk 



Zk~ / ^ \ 2k Zk" Lk" 

— TT^fc \^Kj^| + Lj) j + 2K^g(fc + — gCTj-fe + — (p + P) CTjfe + ^/ {uJk)a + o-ft>a) , (58) 



2k- 



2k^ 



e 



4f , 



e 



= T^kj) - S[kj) + DikSj) + \^K--j £kj + K£kj + —akj + [2K + - j J^kj 

-2'^ik^j) - % (^2^''') " ^'') ~ ^k)) + ■^{0 " i^k)a + <^k)a) - 2% {Aj) - -^Aj) 

I ^ \ ^ 3 K^ _ 

-£{j ° (Wfe)a + ^k)a) + £{k '^'"^j)a [Kb - 2Kbj + £ (^j "■^'Hk)aAb + -j^S l^j"^ CFk)aH-b + -^''^'{kj) 

~2 ~2 ^ ~2 ~2 

-y^(fe^j> + yQ^fej - y (tt - p) S=fej - —qcJki 

~2 r,~2 ~2 



(59) 



= ■H(fcj) + £ah{kD°'^j) + KT-Lkj + QUkj - -^T-L{kKj) - -jS{ji^k) - {^a(k + ^<^a{k) 



+3£(^jQk) ir£ab(kD°''Kj) - K {q{j(^k} + '^{Pk)) ir^(k {^j)a1b + ^j}aT^b) 



e 



4£ 



e 



^(kj) - D^kSj) + [K + -\£kj+ KSkj + —Gkj +\2K--\Fkj- ^-R^k^j) 



^ / ^ \ 3 ~2 ~2 ~2 



~2~ —2,--. _ 



1 



e 



26, 



= H(fc,> + e„6(fcD"J-;> - 2^0-Wfc> + I ^ + 1 Wfcj + —Hkj - 2Haik^^f 



{^k)a — <^k)a) + l^^U^k) — '^'^{k^j) — £ {k ° £j)a [Kb — Lb 
--j^ik^j) + ^(fe "'^^j)^ i^Kb + ift) + -e^fe "'^(Tj)a£b - £{k °'^^j)aA 

o 



49 



g ffc - efe "*?^caJ?,^ + ^£k ^'n^Ab + — - —D^q 



2lf ^ 2lf ^ 2^2 9 , 



-Qk 



2k' , , _a - \ ~„ 2K^g 2k' _ 



+OT-Lk - -£uJk - e 



ab 



k 



-AaSb — -SaKi — La£b 



ab 



DaQb 



ab 



+ D'^Ska - ^DkS + eSk + Y^fe - SkaA'^ - mkaui" +[k- -e j Sk + Tka - 2l» 

-y Dfc (p + TT - p) + ^ (^fc - 2ife) + — ( ^ + -3 ) * - — ^fca^ 



1 _ 25^ 29 

~2 .~2 ~2 ~2 ~1f\ ~2 

+y-D(fe7fj> ^7f(fc^j) + y go-jfc + y (tt - p) ct^tc g^'^J'^ ~ y^C^" ('^■'■>« + ^J>«) 

3 ^ ^ ^ 

~2^{fik) + ^{k^j) - '^'^aik'^j) + ^(j" {'^k)a - ^k)a) + ^{k " ^j)a {Kb + Lb) 

3 ^ ^ 

~^{k " ^j>a i^^b - Lb) - -Sf^^. " aj)a^b + S{k " ^j)aAb , 



'^'{kj) + £ab{kD°'£j) - -j^{kKj) - KHkj + QT-Lkj - {l^a{k + 3(Ta(fc) "H^" + -£ 



U^k) 



-^(k {^i)a - '^^j)a) A + -^£{k "'''^j}a£b " £(fc ''''(^j)a^b " £(fc "''•?\7)a (^6 " '^Lb) 



2 2 



+2 



= ij'^-Wfe - HabCT"'' + UabS"^ - 3^aW" - 35aS" , 

= L)°'£ak — L)"'£ak — '^^k^ ~ '^^k — -^^k + i'Hka'^'^ — £fe '^^'^i-ac'^b" 

{3Qka + ^ka) + (3a;feo + a ka) - 3:Hka^^°' + £fe "'^T-l-ac^b " -^L)°-T^ak 

K _ ^ ^ 2/v O — '2k — 2'-Y7 ^ ~n — 2 — ^ — 

-y -Dfc (P - TT + p) ^^fc ^TTfc + K q Uka + y « + K TT Wfea + —IT Gka 

= D'^Hak + \ek ""'Da^b - ^Hk - Y^k - UkaL" - ^Tka^" - {Uka ' 22^) ^» 

1 

— i^ka — 3^fco) — ej, {Sac — ^ac) b " 2^*: ^aLb + yEfc DaQb 

K „U — 2/t 2/t , — K — K nU — Q 

-Ek TTaLb + -^qUJk + -^{P + P)i^k- y TTfeaW - y £fe 7r„ CTfcc , 

= D'^'Hak — -^Ek '^^DaSb — y^fe ~ HkaL'^ + 2J^fcaW" + ^£fca — 3£fca j 3"'^'^ 
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IV. 3+1 GRAVITATIONAL DYNAMICS ON THE BRANE 

In this section we consider distributional energy-momentum tensor sources on the brane, in addition to the regular 
energy- momentum tensor Tat- Such a distributional source comes together with a discontinuity in the extrinsic 
curvature, as related by the Lanczos equation. 



A. The Lanczos equation 

The extrinsic curvature of the brane is '^^^Kab ~ '^{ci^d) — 9(c9d.)^ '^^t,, equal to the symmetrized version of the last 

four terms of Eq. (I6bp . Replacing Kab by the expression and specializing to the brane cf. Eqs. (ITT|) . the extrinsic 
curvature is expressed as 

'-'^^Kab = KUaUb - 2U(^aKb) + ^hab+ ^ab ■ (74) 

As we approach the brane from left or right, the limiting values of the extrinsic curvature could be different, according 
to the embedding and 5d metric in the two regions. Therefore we introduce averages and differences of the extrinsic 
curvature. 

The Lanczos equation [46|, relates the jump of the extrinsic curvature across the brane to the distributional 
matter layer: 

A^^'^Kab^-K^ [rab-^gab] ■ (75) 

The u°-u^, u°'h'^, trace and trace-free parts of the h'^h^ projections give: 

~2 

AK = ^(A-2p-3p) , (76) 

AKa = K^a , (77) 

AG = -K^i\ + p) , (78) 

ACTq6 = -K^TTab ■ (79) 

The Lanczos equation is necessary in order to derive the gravitational dynamics on the brane, given by a scalar (the 
twice-contracted Gauss), a 4d vectorial (the Codazzi) and a 4d tensorial (the effective Einstein) equations Q. The 
latter has been first derived in [48| , later generalized to include bulk matter and asymmetric embedding contributions 
(Eq. (1) in i3|]). 



i 



B. The 3+1 decomposition of the source terms of the effective Einstein equation 

We give the 3+1 covariant decomposition of the source terms of the effective Einstein equation. This equation is 



Gab = - ( ^ - '^-^ ) 9<^b + Ab + H^Sab - l^-^Kb) + {Lab) + {Vab) ■ (80) 



The sources are: the stress-energy tensor Tab representing standard model matter [decomposed in Eq. (flSl) ]: the source 
term Sab quadratic in Tab (dominant at high energies), {Vab) the pull-back to the brane of the bulk matter; {Lab) a 
source term originating in the asymmetry of the embedding and ^^■'fab the contribution of the electric part (relative 
to the vector n") of the 5d Weyl tensor. We have defined the 4d coupling constant and the brane cosmological 
constant A as 

Gk^ = k'^X , (81) 
2A = k^A+Za) . (82) 
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The quadratic source term is decomposed as 



Sab = ^ (2p^ - SlTcdTT'''^) UaUb + ^ {2p'^ + 4pp - Aq^q" + T^cd'^'"^) hab 



+ ^q(aqb) + -^q(aUb) - -q^T^daUb) 



P + ip 1 
—^T^ab - -^T^cia-^b) 



(83) 



Some of the numerical coefficients here are corrected with respect to the corresponding expression (7) given in Ref. 



The electric part of the 5d Weyl tensor expressed in terms of gravito-electro-magnetic quantities defined in Section 
HTCl is 



The asymmetry source term is decomposed as 

{Lab) = : 



UaUb 



+ - 



-U(a 
1 



b)c 



2 {Sh)c) 



e 



+ {Ka){Kb) + (--K) (dab) - {da') {^bc) 



(84) 



hab 



(85) 



As the induced metric is continuous, the average of the trace L ~ g Lab is the trace of the average: 



{L) = {dcd) (r') 2 (Kb) (k") ~ I (§)' + 2 (e) (^K) 



(86) 



For a symmetric embedding {^*'>Kab) = 0, therefore cf. Eq. ^ (^ej = (^Kj = 0=(^KcJ = {^cd), so that (Lab) 
Finally, 

6Vab 



= 3{p+p) UaUb + 8q{aUb) + (p + p) hab + All ab 



(87) 



C. Gravitational dynamics on the brane: generic embedding 



In order to obtain the evolution and constraint equations on the brane, we select a subset of the Ricci and Bianchi 
equations given in subsections IIII Al and IIII C\ by combining them in such a way, that the off-brane derivatives of the 
kinematical and gravito-electro-magnetic quantities drop out. Additionally, the equations of this subset contain only 
quantities appearing in the effective Einstein equation and in the 4d theory. 

First we express Ha, £a, ^ab and "Hab from ((37)) . ((50)) . p3)) and p9|) respectively and we employ the definitions p4| . 
()15|) in order to introduce Eab and Hab in place of £ab and T-iab- Inserting these into the system of equations given 
by the following equations: dHD, O, O, (EH)-®, (TO, (ED, (EHD, (ES), 60]), O-dnil), 0, GD) and 

evaluated at the brane, we obtain a system of equations to be referred as the brane Eqs. These equations are either 
evolution or constraint equations on the brane and for a generic asymmetric embedding are presented in Appendix 
[Cl The evolutions refer to the quantities 8, Ka, O, Wq, Oab-, ^a-, Eab-, Hab- 



D. Gravitational evolution and constraint equations on a symmetrically embedded brane 



In this subsection we restrict ourselves to symmetrically embedded branes. The Z2-symmetric embedding arises 
when there is a perfect symmetry between the 5d space-time regions on the two sides of the brane. In this case the 
extrinsic curvatures on the two sides of the brane are opposite. This is due to the fact, that the normal vectors to the 
brane on its two sides are and — respectively. Therefore A^'^^Kab = 2''^^ifab and Kab = 0. 
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We present a system of evolution and constrain equations, which hold on the brane and contain no off-brane 
derivatives .We obtain these equations by specifying Eqs. (|Cip - (|C14p for a symmetricahy embedded brane; then 
replacing Q, Ka, K and aab with the corresponding matter variables as given by the projections ([75|) - d75|) of the 
Lanczos equation, again specified for a symmetrically embedded brane. Finally, we employ the definitions ()81|) and 
(15^ . whenever possible. 

For improved clarity we group all general relativistic contributions on the left hand side of the equations, keeping 
the brane-world contributions on the right hand side: 



p + D^qa + ip + p)e + 2q''Aa + TTabCr 



ab 



46 



q{a) + DaP + D -Kab+ ^<7a + {p + p) Aa + TlabA - U abq + Cf abq 



-2q 



-2n„ 



(89) 



= 



TT TT, 



+ ^'^'Daqb + q'^D'TTab - ^q^'DaP 



-2Abqa7r^' + l^ah^"' +{P+P) '^abTT^' + ^caTT, + ^qaQ'' " 'T.^q'^q'' 



29 

-4 



- Y (P - ^ + - - ^© iP + P} + 2K'g + ypq - ^'^'^^^ " ^^a&CT°^ , 

£(k) + 3®'^fe ^ Z^^^ 3"^*^ ~ D°-^ka — S-kaA"^ — (Wfea — d ka) 



K 



29 



-'T{/ca)9" - iP+P)D TTkb - q^CTckT^^" + - {p + p) Dkp - — [p + p) qk 



1 ^. 2 y. 

'2q''Dkqa + q^Daqk + ^QkD^'qa ~ 2q°-Ai^kqa) + (^baT^k ^ ^qk(T abT^"-" 



K 

y 



1 K 

-2i'^fe) + -T^fc (p + 3^ - 3p) + y (2A - p - 3p) - 2Xgfc 



y(7 (2^fe + Xfe) + 27rfe„> - 2 (tt - p) A. + ^TTfeaTr" 



(90) 



(91) 



9 - D'^Aa + y - ^"^a - 2^,0.'^ + (Tafcfr'^'' + y (p + 3p) - A 

~4 ~4 ~2 

^ + yOa-^p(2p + 3p)-y(p + i+p) , 

n . 29 . . ^ ^ 
"■(ab) - -0(0 A) + ^CTab - ^(aA) + ^(a^h) + c{a<^b) + -^ah ;^7ra6 



(92) 



K 

"3 



1 29 



(93) 
(94) 



D'aab - ^DaQ + £„ "'^'I^cWfc + 2e^ ""A^u^k + n'qa 



-^pqa + y VTo 



(95) 



£'(c'^fc) + SabikD^a^f + 2A(,a;fe) + Hafc = , 

D'^UJa - AaLo'' = , 



(96) 
(97) 
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K 



e 



Hkj) 



a ' ^ j)a ) 



2^{kj) - 2^{k^j) + 2^0" i^k)a + <^k)a) + -^^kj - ^{k^j) + "^^^fcj 



K 

24 



{p + 3p) 7T(fej) + 67r//7Tfe)a + (p + ^p)' TTkj " 2g(fc£)j) (/3 - 3p) - 2pD(kqj) + 3Tr,j'' Dk)qa 



+6q^kDKj)l, + SqaD^kTrj)" -2p{p + p) akj + l^q^k^j) + 2 (p + 3p) g^fc^^) + OTr^/TTfc),, 
+ ^ (P + 3p) TTfcj + (p + 3p) TT^fe" + (7j)„) + 'iq(k<^j)hq^ + 37r^j"a;fe)c7r„'= + SajkQ^Qa 



+- 



TT{kj) + Dl^kQj) " '^'^{klj} + 4a(fc-4j) + {p + p)ajk + —TTjk + TT^/ {uJk)a + (^k)a) 



(98) 



2 

H{kj) + £ab{kD'^E^'' + OHkj - 3cra(fci?j)" - U:a{kHj^°' - 2e^^ "■''Ej^a^b ^ {sahikD^'TTj^ + -ioj^kQj) + £(fe °''''^j)aQbj 



ia_ 6 



'=''7r,>e£>d (p + 3p) - (p + 3p) s^bikD^n^ 



Se^bikD^TTj^^w"'' - 3eab{kTTj)cDV + Gpu^kqj) + "^pSik "''^j)aqb + Se^b{kQj)D''q'' + 3e^b{kQ''D''qj) 



'\a—cb 



■)0 -,6 



j) 



+3£,fe "''aj^bTTa'^qc - QTT,j°'UJk)qa 



+ y {£ab{kD°-Tij) + Sq(_jUJk) + e^^ "'''crj)aqb) 



(99) 



D°'Eak — ^Hka'^"' + £fe "'^Hac^b 



2 2 



-pDkP + TTafe-D" (p + 3p) 



I'^l ©'^l '^K 

-D°'£ak — '^E>k£ — —£k + 2 (3w;sa + <''fea)£° + ^ 

+ (p + 3p) D^TTak + 37r/£»°7refc - 47r»''A7r„6 + 37v kbD^n^ ^ + 2q"Dkqa - 3qaD"qk 
-3qkD"qa + 96^ "Sc^aTTd" - 37r/gVfe6 - ^pQqk + '^pc^kbq^ + ^P^k ^Uc^d + '^^'^k'^qa 



~2 
K 



1 „ ~ 2 „ 

i'^TTafe - -£>fe (p - TT + p) + -@qk - q" {Su^ka + (^ka) 



(100) 



D'^Hak + 3£toW'^ - £fc '^''£;„eab " + ^ [£fc ""-Dagft + 2 (p + p) Wfc - TTfe^Wc - £fc ""TTaca;, 

Sq^oJcqk - 2£fc ""^qbDaP 

-2p£fc ''''iP^gft - 3£„6feg'=L'''7r, " - 3£fc "''w^^Dbqc - 4p (p + p) Wfc + 35ag'^Wfe - (p + 3p) tt^^Wc 
-3£fe "V^feg'-ge + 37r„67r«*a;fc - 37re„7rfc'=w» - (p + 3p) «V„c^6' - 3efc "^daTr, '^a^, ' 



5?fc„a;" - £fe 56 - 2 (p + p) wfe + £jt ""S^,, 



(101) 
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For vanishing 5d matter these equations reduce to the corrected form of the Eqs. (26)- (29) and Appendix A of Ref. 

m- 

Eqs. ()88|) and ()89|) express the interchange of energy density and energy current between the brane and the 5d 
space-time (due to the nonvanishing right hand sides). In the absence of the 5d sources, these equations become 
evolution equations for the brane matter alone. Similar relations for the effective nonlocal energy density and effective 
nonlocal energy current are given by Eqs. (j90| and (|9T|) . 

Even for the chosen Z2-symmetric embedding and with the simplifying assumption Tab = (thus = ]) ~q = qa = 
TT = TTa = TTab = p) the abovc cquations do not close on the brane, due to the absence of evolution equations for Sab, 
TTab, and p. Assumptions fixing Wab (typically by kinetic considerations employing the Boltzmann equation) and p (by 
choice of a continuity equation) are required already in general relativity for closing the system, however on the brane 
(even with empty 5d space-time) an additional assumption for Sab is equally required jTlj . From these considerations 
it is immediate to establish the closure in the special case Sab = = iTab, provided the equation of state is known. 

As the closure is difficult to achieve even in the simple 5d vacuum case, in order to tackle realistic problems we 
need to consider the ensemble of all dynamical and constraint equations given in the preceding section. 



E. Closure conditions 



The general relativistic 3-1-1 covariant formalism contains 10 gravito-electro-magnetic variables with 10 evolution 
equations (besides there are also 12 kinematic variables with 9 evolution equations and 15 constraints altogether). 
The 3+1+1 brane-world formalism developed in this paper contains 35 gravito-electro-magnetic variables with 35 
evolution equations (beside 35 kinematical variables with 28 evolution equations and 77 constraints altogether). 

The subset of equations on the brane, given in subsection IIV Dl has 10+9 gravito-electro-magnetic variables with 
only 10+4 evolution equations (there are also 12 kinematical variables with 9 evolution equations and 15 constraints 
altogether.) The 9 new gravito-electro-magnetic variables are the quantities appearing on the right hand side of Eq. 
(|84p : among them the last term Sab, representing 5 independent variables, has no evolution equation. It has been 
known that the system of equation is closed by the condition Sab = [Tlj . 

In this subsection we want to explore the extra information we have in the complete system of evolution equations 
derived. In particular, Eq. (j6ip contains the desired temporal evolution, however it also contains terms not appearing 
in subsection IIV Dl Remembering, that on the brane Qab = we could impose 




+ pikAj^ - ''%)a {Kb - 2kb) - ^''nu)aAb + "''?fc)aHb , (102) 

such that Eq. (|6ip becomes 

= Si^k]) - D(k£j) + ^£kj + yCTfcj - 2f (fc^j) + {ujk)a + eTk)a) + Mkj , (103) 
with the 5d matter contributions 

~2 ~2 ~2 ~2~ ~2r\ 

o~2 ~1 , , ~2 

+ y'^'O" [^^k) + Lk)j + yT^o" {^k)a + cfk)a) ■ (104) 

A particular solution of Eq. (|102p would be 

JTj^^- = -Hfej = Sj ^ Uj ^K = Aj=Q. (105) 

None of the quantities above appear in the brane equations presented in subsection IIVDI thus those equations are 
not altered by the choice ()105p . and the system becomes closed by Eq. ()103p .'^ 



^ The no-go theorem for closure given in Ref. |32ll does not apply here, as it was derived for perturbations of 5d Anti de Sitter spacetimes. 
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The 5d Schwarzschild-Anti de Sitter space-time containing a Friedmann brane emerges as special case of the space- 
times obeying Eqs. (|103l) and (jlOSp . as they have Mat = — K and 

£a = £a = = La = Ka = Ka = UJa = i^a = Aa = Aa , 

= Sab = J-ab = "Hab = £ab — ^Hab = ab = ^ab ■ (106) 

At the end of this section we emphasize, that Eq. (I103P closing the system of brane equations should allow for far 
more solutions than the trivial one for Sab- We will construct a specific example in Section IVll 

V. COSMOLOGY 

In this section we consider generic embeddings. By employing the definitions (|8T|) and ((82|) . also the conditions 
(jlip . arising from the existence of the brane, we can derive average and difference equations on the brane. We give 
here but the most relevant such dynamical equations. The rest of the equations is straightforward to derive, although 
it may be lengthy. 

The average taken on the two sides of the brane of Eq. ()25|) reduces to the generalized brane Friedmann equation. 
We also rewr ite (e^'j = (e'j + (ag) /4 and use a similar relation for dab- We take AQ from the Lanczos equations 
(US])- 111]). Eq. ^ then becomes 



| + ^- A-«.V(l + ^) - I'^aba'^' + ^aUJ^ = -{£)- ^nabn^' + ^ - I {dab) {d'-') + y (p + p - ^) . (107) 



The generalized brane Raychaudhuri equation is obtained from the average of Eq. (j27p . by employing the same 
sequence of simplifications, as for the Friedmann equation. We obtain: 

e + ^ + aabcj"'' - iujaoj" - D'^Aa - AM„ + ^ (p + 3p) - A 



{£) - (2p + 3p) + yga?"^ - (e) (k) + (ka) (^k^) - ^ (p + n + p) . (io8) 



Finally we give the generalized brane energy-balance equation from the jump across the brane of the Eq. (|4ip by 
the same procedure: 

p+e{p + p) + Daq" + 2Aaq'' + Tiab<y^^ = - Aq . (109) 

These equations hold for arbitrary branes. Again, general relativistic contributions are on the left hand sides, brane- 
world contributions on the right hand sides. 

A. Friedmann brane with perfect fluid 

In this case the conditions uja — — Uab = Auab hold, arising from the particular geometry and matter source 
chosen, also TZ = Qk/a^, where a is the scale factor. Moreover 

Q/3 — H = d/a, where H is the Hubble parameter. The Friedmann, Raychaudhuri and energy-balance equations 
become: 

3 ( H 



2 



{dab){d'"') + Y(p + p~^) ■ (110) 



2 

^4 _4 



3 (i? + + ^ (p + 3p) - A = {£) -'\^i^P + 3p) + ^qaq" 



e)(^K) + (^ka){K^)-^{p + n+p} , (111) 
p + 3H{p + p)^-Aq. (112) 
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For symmetric embedding the equations simplify by = = (ffab)- For generic asymmetric embedding Eq. (jllOp 
can be rewritten as 

where (L) is given by Eq. ([55]) and U is defined by 

^'U=^ + \ (e) (i?) - i (i^,) (X^) - i (a„,> (d'^") - (£) + I (p + p) . (114) 

The quantity U is nothing but the effective energy density introduced in Ref. encompassing the trace- free parts 
of the Weyl, embedding and 5d matter sources in the effective Einstein equation, given here in terms of the 3+1+1 
kinematic, gravito-electro-magnetic and 5d matter variables. 



B. Anisotropic brane- worlds 



Full brane-world solutions with homogeneous, but anisotropic 5d space-times are also known. 

In Ref. pgj a vacuum 5d static and anisotropic space-time with cosmological constant admitting a moving Bianchi 
I brane was analyzed. The Z2-symmetric junction conditions could be obeyed only by anisotropic stresses on the 
brane, hence the brane cannot support a perfect fluid. Isotropy of the brane fluid could be achieved only for isotropic 
5d space-time and brane. 

This setup was generalized in Ref. [50| by allowing for a non-static 5d space-time. Assuming separability of the 
metric components, new 5d solutions combining the 4d Kasner solution with the static 5d solutions of Ref. 49] were 
obtained. 

For the reader's convenience we give in Appendix |D] the list of kinematical, gravito-electro-magnetic and matter 
variables for the 5-dimensional models presented in Ref. [49| . Working out the respective quantities for other metrics 
would be a similar straightforward exercise. 



VI. STATIONARY VACUUM SPACE-TIMES WITH LOCAL ROTATIONAL SYMMETRY 

In this section we discuss an application of our formalism, by assuming vacuum in both 4d and 5d, but keeping the 
respective cosmological constants. The embedding of the brane is symmetric. Then the effective Einstein equation 
reduces to 

Gab = Agafc - f [uaUb + ^^^ab^ + iSf^aUb) - Sab ■ (115) 

We are interested in stationary space-times, therefore / = for any scalar field /. The stationarity implies a 
singled-out temporal Killing vector, therefore the 3+1+1 decomposition of the gravitational dynamics developed in 
this paper turns particularly useful for the study of gravitational dynamics on the brane (which defines the other 
singled-out direction). 

We also specialize our treatment to space-times with local rotational symmetry (LRS) on the brane. Such a 
symmetry singles out an additional space-like unit vector field e", in the sense that there is a unique preferred spatial 
direction at each point that assigns the local axis of symmetry. Once such a special vector field is chosen, a further 
decomposition of the spatial quantities would lead to a generic 2+1+1+1 formalism. For this, the metric hab should 
be further decomposed as 

hab = SaCb + qab , (116) 

where qab is the induced 2-metric on the surface perpendicular to both e'^ and u"", and lying in the brane. 

In what follows we will see that these symmetries assure that the structure of the space-time can be described solely 
in terms of scalars, thus no vectors and tensorial quantities are needed. This is a powerful feature of the formalism. 
Furthermore, the symmetries assure, that all scalars / depend only on the coordinate parametrizing the integral 
curves of the rotation axis field e°. We denote this coordinate derivative as /* = eJ^Daf (a spatial covariant derivative 
along these integral curves). 
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A. Independent kinematic quantities related to the vector field e" 

1. Decomposition 

For the purposes of the present application it is enough to give here the decomposition of the covariant brane 
derivative of the vector field in terms of kinematical quantities and extrinsic curvature type quantities analogous 
to the ones appearing in the decomposition of the vector fields u"- and n"- . We keep the familiar notations, with the 
remark that the quantities belonging to the decomposition of V a&b will carry a distinguishing~mark. We thus have 

VaCf, = LUaUi, - LaUb - UaK^ + Ba^fc + ^Qab + ^ ab + ^^ab , (117) 

with 

La = u'^ha^VcCd , Kb ^ u'^h^'^V cCd , A = e^^g^ ''VcCd , 

Wab = q^a''lb]'^^ced , '^ab^ qa'^qb'^^ced-^qab ■ (118) 

We remark that L and La are not independent from the previously introduced sets of variables, they can be expressed 
in term of projections of the kinematic quantities related to u": 

L = -e''Aa , La^-e'' (^^had + CTad+i^ad^ ■ (119) 

By contrast, the quantities Q, Kb, Ab, Uab and ujab are independent of the rest of the kinematical and extrinsic 
curvature- type variables. Similarly to and Wq, we can also define a rotation vector 

l^a = '^Sabc^^" ■ (120) 

2. LRS symmetry 

The preferred spacelike unit vector field Ca satisfies m^Cq — 0, e°-ea — 1- We employ here various results following 
from the LRS symmetry, following Ref. (sij . The symmetry and normalization implies: 

e^Daeb = , e(b) = , (121) 

i.e. Co is geodesic with respect to the connection compatible with hab and is Fermi propagated along the world line 
of a brane observer. The above equations and normalization further imply that Eq. ()120p can be rewritten into the 
standard form 

= Sa "^gib ^g,] "^eed = SabcD^e^ . (122) 

Due to LRS, all spacelike vector fields must be proportional to e°, thus 

= Ae°- , = WCa , UJa = ^^a ■ 

r = i^e'' , i?°e = e*e„ , Da£ = £*ea , (123) 
The vector field e° and the induced metric qab define a unique spacelike tracefree symmetric tensor field eab as 

eab = BaCb - ^ , (124) 

satisfying: 

u'^Eab = , e^Cab = Cfc , e'\=0 , iCacE^'b = C^ab + hab , 

3 30 

eafeC"^ = - , D''eab = -Y^a , e^ab) = . (125) 
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Again, due to LRS all 3d tracefree symmetric tensor fields are proportional to eab- 

2(j 2E 2H ^ 2£ 

In Eqs. (|123p and (|126p we have introduced suitably normalized scalars A, w, w, cr, £^ , £, E and iJ, replacing the 
vectorial and tensorial variables. 

From LRS, by use of Eqs. (fT2T|) . (fT23l) . (fTM]) and (fT25l) also follows that 

L = -A , La = -(j + ^ea , (127) 

i^a = la = , 5afc=0 . (128) 

Thus there arc only two kinematic quantities related to e"' left, which are both non-trivial and independent from those 
introduced the Sec. [Ill These are 8 and uj. 



B. LRS class I type conditions 

The general relativistic classification of the LRS models presented in [51] is recovered for £q = = £ab and £ <Q. 
For brane-worlds, when the above conditions do not hold even under the simplifying assumptions of this chapter, this 
classification should be refined, however for the application we are interested in, we shall still assume the conditions 
cj = 9 = cr = characterizing the LRS class I of general relativity. From these conditions, Eqs. (^51) . (|122l) . (|123p 
and EaijC^e^ = we also find £^ = 0. Therefore we have verified, that Eq. (|103p . which closes the system of brane 
equations for the considered stationary vacuum space-times containing LRS class I type branes, trivially holds. 



1. Dynamics 

The evolution of the single kinematic quantity characterizing e° can be inferred from the Ricci identity for e°: 

(129) 



~ _ 2E 2£ £ 2A 



Other nontrivial brane Eqs. are (l9T|), ((921), dMl, (EH), (El), (flOOl) - and (fTOTj) . Under the assumptions of this section 
they simplify to: 



£* + AA£ + 2\/3 



A* + [Q + A] A + 2u^ + K 



A* + \ A 



A-u? - ^fZE 



, 
V3f 



ij + 2A - e 



= , 



e-A w = , 



^ 36 £* ZQ£ £* 

E" + —E-3ljH = — + — p 

2 2 4 2^3 



2H* + mH + &Euj ^ 



Auj£ 



— uj£ 



Eliminating A* from Eqs. (|13ip and (|132p we obtain the algebraic equation 

= ^+3c.2+A-fV3i? + £ + ^ 



(130) 
(131) 
(132) 

(133) 
(134) 

(135) 

(136) 

(137) 
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Also for any solution of the system (|129p - (|134l) . Eqs. (|135|) and (|136p are identically obeyed. This can be seen as 
follows. By taking the ★-derivative of Eq. (|137p and employing Eqs. (|129p - (|134p we obtain Eq. (|135p . Similarly, the 
★-derivative of Eq. (fTM]) . combined with Eqs. (fT29l) - (fT3tl) gives (fnej) . 

As we have 6 independent (2 algebraic and 4 first order differential) equations left for the 7 variables (£, 8, 0, A, 
Lo, E, H), we need to impose an additional ansatz, chosen as 

- IF 

£ = -^. (138) 
V3 

This condition will considerably simplify the algebraic equation Eq. (I137p . 

2. Discussion 

The algebraic Eqs. (|133p and (|137p give H and E in terms of the rest of variables. By Eq. (|138p the system 
(IT ^ - p?T|) and (UMl) reduces to 

~ ~ 

+ — - QA-2uj^ ^ , (139) 

£* + 2£Q = , (140) 
A* + {q + A + 2uj^ + K = -£ , (141) 

uj* A^uj = . (142) 

From the newly arised two algebraic Eqs. (I140p and (|142l) we express 

e = {\n£-^'^y , (143) 

A = (inwf-i/^y (144) 

In terms of the auxiliary variables 

X = \n£-^'^ , (145) 

y = Incjf^i/^ , (146) 

the remaining Eqs. p39p and (|14ip become 

X** + - x*y* ^ 2ey-^ , (147) 

y** + {y*f + x*y* = - (2e«-^ e'^^ A) . (148) 

They form a coupled second order system, which would eventually give w and £ in full generality. The solution of 
this system is however not immediate, therefore in what follows we will employ a metric ansatz in order to find a 
particular solution. 

C. Taub-NUT-(A)dS solution with tidal charge 

We take the following metric ansatz, compatible with the chosen symmetries: 

ds^ = -l^ {dt + uj^d^f + ^^dr^+g (r) idO'' + flld^^) , (149) 
g{r) f{r) 



where 



sin^ , k = 1 

r^k = <( 1 , k = , (150) 
sinh6' , k — —1 
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and is another function of 9. The axis of LRS is given by 



Choosing the 1-form as 

fl/2 

u= -^{dt + 2uji,dif) , (152) 
and employing £^ = 0, Eqs. (|126p . (|138p . (|151l) . the electric part of the 5d Weyl tensor is found as 

= (4)£:-^ = - i*)£\ ^ _ (4)f ^ „f , (4)£:*^ ^ _2£:wk • (153) 

Both from £a £^ = and from Hat being proportional to Cab we find 
Equivalently, by an integration: 



n^'^^-2l (155) 



where Z is constant. A second integration gives 



2/cos6' + L, k = l 
Wk(6')=<( -210 + L , fc = , (156) 
2/cosh0 + L, fc = -1 

where L is another integration constant. Locally, this constant can be absorbed in a new time variable t + Lip — >■ t, 
which translates to the choice L = 0. 

Direct computation, employing Eq. (jl56p shows that the metric ansatz is compatible with the chosen LRS class I 
conditions 

Q^a^a = £^'=0 . (157) 
The nontrivial kinematic and Weyl quantities, by employing Eqs. (|149p and (|156p are 



e = 


/V2 

g3/2 ^J. ' 


dr 




55/4 dr \ 


53/4 rf^y 


2f 




2 rf/dg 




3g dr"^ 




2^; 




1 d 


7! ~ 


6(7 rfr^ 


dr ' 


2H 

71 " 







5^ 5 



4P/ 2k 
H ^ H 

f/^dg\ APf k 



(158) 

These quantities are constrained by Eqs. (|133p . (|137p . (I138p . and (|139p - (ll42p . From here by straightforward algebra 
we find two independent equations for the metric functions / (r) and g (r): 

-4 ^ 2k- 4Ag . (160) 
dr^ 
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By multiplying the first equation with dg/dr and integrating, we get 



2 



41' = Cig , (161) 



with Ci > an integration constant. A further integration then gives 

C 4/^ 
9{r) = ^{r + C2f + -^^ , (162) 

with C2 another integration constant. The constant C2 can be absorbed into r by a redefinition of its origin, therefore, 
without restricting generality we choose C2 = 0. Also the constant Ci disappears in the following rescaling of the 

coordinates and parameters: (2^/Cy^ Cj/^r/2, 4P/C1) ^ (t, r, P). Formally this is the same as choosing Ci = 4. 

Eq. (|160|) then gives the other metric function as 

/ (r) = C4 + C:ir + {k - 2fA) - , (163) 

o 

with C3 and C4 emerging as integration constants. With the reparametrizations C3 —2m and C4 = g — fcZ^ + A^^ 
we find the brane solution given by Eqs. (|149p . (|150p and 



f{r) = k {r^ - f) - 2mr + g - A + 2Pr' - l"^^ 
g{r) ^ r'^ + l\ 



2Zcos6i , k^l 

LOk{e) = { -2ie , fc = o . (164) 

-2;cosh6', k = ~l 

This is quite similar to the charged- Taub- NUT- (A)dS solution of a general relativistic Einstein-Maxwell system with 
mass m, electric charge Q, NUT charge Z, and cosmological constant A, however the constant q replacing Q' is not 
restricted to positive values. A glance at the effective Einstein equation (IllSp shows that this constant could possibly 
arise only from the electric part of the 5d Weyl tensor, Eq. (jl53l) . Indeed from the fourth Eq. (I158P we get 



£ = ^ . (165) 

(/2 + r2)2 

As q originates in the Weyl sector of the higher dimensional space-time, the derived solution has the interpretation 
of a Tauh-NUT-(A)dS brane with tidal charge. 



VII. CONCLUDING REMARKS 



In this paper we have developed a generic 3+1+1 covariant formalism for characterizing 5d gravitational dynamics 
on and outside a brane. The singled-out directions are the off-brane and temporal directions, thus the 3-spaces are 
constant time sections of the brane. Generalizing previous approaches, like 3+1+1 with the extra requirement of double 
foliability [Qj, [10!|, 3+1 in general relativity [l2|-[l^ and in brane-worlds [ll], finally 2+1+1 in general relativity jl6j . 
p~7l |. we presented gravitational evolution and constraint equations in terms of kinematic, gravito-electro- magnetic 
and matter variables, defined as scalars, 3-vectors or 3-tensors. The number of variables being higher than for other 
lower dimensional formalisms, the 5d matter and especially the 5d Weyl tensor leads to a multitude of projections 
without counterpart in the mentioned approaches. Only the kinematical set of variables is similar. We have compared 
and checked our results with those presented in a recent work by Nzioki, Carloni, Goswami and Dunsby [isj on the 
2+1+1 decomposition of f{R) gravity and we give the correspondence between our and the notations of Ref. [l^ in 
Table 1. 

Our generic formalism contains the complete set of dynamical equations in the 5d space-time. All equations, with 
the exception of those in IV to VI are independent of the particular form of the dynamics on the brane. As such, they 
can be employed both to discuss DGP / induced gravity type jisj branes (a project deferred for future work) and 
one-brane Randall-Sundrum type branes. For the latter, we have given both the full set of equations on the brane, 
and those containing off-brane evolutions. 
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TABLE I: Comparison of the notations for the kinematic quantities with Ref. [l^ . 
Ua^ea /C->ie + E K ^ -A 

Aa Aa Ka ^ Sa + Eab^^ R'a da 

Aa-^aa e^e-|E Q^(tl 

La 'Sa ~ Eabfi'' i^ab f^Eab i^ab ^ ^£ab 



hab ^ Nab ab ^ Sab ab ^ C 



ah 



The brane equations are more general than previously published results by the inclusion of arbitrary 5d sources. 
Although we have recovered the known fact that in the generic case this system does not close on the brane (excepting 
the trivial case Eab = 0), by deriving the complementary system of equations of the 5d dynamics we could establish a 
more generic condition for closure. This is given by either of Eqs. (|102p . (jlOSp which carry much richer possibilities 
than the previously known £ab = case. 

The initial value problem in general relativity is usually discussed in ADM-like variables (including modifications 
enhancing stability, like the use of variables with factored-out conformal factors j54(]), therefore a similar treatment 
would be possible to develop in the framework of the Hamiltonian approach presented in A 3+1 covariant 

approach for discussing the evolution of cosmic microwave background anisotropics in the cold dark matter model 
was employed in Ref. [24| . The complete set of infinitesimal frame transformations given in the present paper may 
turn useful in the study of perturbations in a 3+1+1 setup. 

We have decomposed the Lanczos equations and all source terms of the effective Einstein equation in terms of the 
3+1+1 covariant variables. The ensemble of these results opens the possibility for applications, with both cosmological 
and other symmetries. 

We have given generalized Friedmann, Raychaudhuri and energy balance equations for a generic brane, and by 
specifying for cosmological symmetries we have established correspondence with previous related work We have 
also established the correspondence with the anisotropic brane- world presented in Ref. (49j . 

We have also employed the 3+1+1 covariant formalism to discuss stationary space-times with local rotation symme- 
try of class I, imposed on the brane. We have shown that such space-times obey the closure condition presented here. 
The symmetries and metric ansatz (jl49p implemented in the 3+1+1 covariant formalism led to a simple decoupled 
system of second order differential equations for the metric functions, the solution of which gave a new exact solution 
of the effective Einstein equation, the tidal charged Taub-NUT-(A)dS brane, given by Eqs. pSOp and (1164 1. 

In the spherically symmetric and rotating cases tidal charged brane solutions were already found [iS^] - 53j . which 
correspond to electrically charged general relativistic Einstein-Maxwell solutions, when a formal identification of the 
electric charge squared with the tidal charge is carried on. Here we also found that replacing the electric charge squared 
in the electrically charged Taub-NUT-(A)dS solution of an Einstein-Maxwell system with a tidal charge originating 
in the Weyl curvature of the 5d space-time leads to a brane solution with the same symmetries. Unless the electric 
charge squared, the tidal charge however is allowed to have both positive and negative values, thus allowing to either 
weaken gravity, or contribute towards its confinement on the brane. 
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Appendix A: Commutation relations 

In this appendix we enlists some useful differential identities, obtained by computing the commutators among the 
derivatives u°''S/ a = D/dr (dot), n^Vo = D/dy (prime), and Da (induced metric compatible covariant derivative) on 
scalars, brane- vectors and symmetric brane-tensors of second rank. 
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For a scalar field (j) the following commutation relations hold: 

n"Va(?i) - w^Vc(<^') = Kcp' + K^+ (k'' - Dacj) , (Al) 

DaCj>' -h,'n''Vb{DiCt>) = -{Ka-La)^ + Aa^' + ^DaCj)+{0Jab + aab)D''ct> , (A2) 

Da^-hJu''^b{Di(l>) = -Aa^+(^Ka + La)(t>' + jDa<j>+{uJab + (Tab)D''(l>, (A3) 

D[aDb](l) = Wab(l>-^^ab(l>' ■ (A4) 

For a 3-vector field V"' the following commutation relations hold: 

-AaVKb + KaV"Ab + KaVAb - AaV"Kb , (A5) 



hjh^ ^u^Va {DiVj) - D, (y^d)) = - (i?c + Lc) ^(d) + ^cT^<d> - ^D,Vd + jKdVc - jA^Vc 



2 



— VcQd -{^hcd + Wed + S=cd^ + Kd {Qca + ^ca) V" 

+ (^hcd + Wed + «Ted) AaV - Ad {Wca + (Tea) V , (A6) 



K'h^ 'n^Wa (A^,) - A (V(d>) = - (ic - ife) V(d> - leV^'d) " -^Vd " ^KdV, + -AdV, 



1 1 ~^ 

- (^ae - 2„e) D'^Vd - -SaV^Kd + ^£dVc - Sdab'Hj'V'' - y^FaV^/led 

+ yK7rd + ( ^/icd + Wed + 0-ed J i^aV^" - Kd (Wea + CTea) V 

- ( |/led + Wed + ^ed J ^aV^" + Ad (Wea + CTea) V , (A7) 



-Dfa-DfejK = Wa6V(c> - WafcV^c) + hc[a£b\dV'^ + ^ela^b] - /lc[a'^^6]d^'^ - ^e[a^6] 

(©^ - 0^) /le[a^6] - ^ (o-e[o " We[a) H] + ^ (S=c[a " We[a) Vi 

--^hc[a {crb]d + Wfcjd) F'' + -^/lc[a (o'hld + Wfcjd) V^'' - K[aVb] 



- {(Tc[a - Wc[a) (Wb]d + Cr;,]d) + (^efa " We[a) (w6]d + 5^6]d) 

T ~2 ~2 ~2 

+ -^hc[aVb] + Y (p - TT + p) /le[aH] + y /ic[a'f5]d^'' + y TTelaH] • (A8) 

For a symmetric trace-free 3-tensor field Tab the following commutation relations hold: 

\c%)'n'^^a(T^ij)) - hoK)'^'^b{T[^^)) = KT[cA) + K%d) + i^^ATed - K-D,T,d 



-2K^,Td)aA^ + 2A^,Td)aK'' + 2A(eI^d>a^" 
c) ^ ' 



-2K^,Td)aA'^-2eab(dT,)'"H\ (A9) 
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ha 'h^.'h^^ 'n'^Vd {DkT,j) - Da (r(,^)) = - {La- Ka) T(bc) - AaT[,^-^ - jDMc + £(^bTc)a 

^ 20 

- {^ad + (fad) D The - 2T^^ ec)di'Ha^ + —ha{bTc)dK 

-ha{bT^)'^£d + -^A^bTc)a ^ha{bTc)dA''- 

-2 {iJad + ^ad) Ki^tT,^ + 2 (Qad + ^ad) 1(6^,) 
+ 2 {uJa{b + O'a(h) Tc)dK'^ - 2 {uJa(b + ^a{b) T'c>d^'' 

26 2k\ 2k^^ 

~-^K(^bTc)a -^ha{bTc) TTd + -g-7r(bTc)a , (AlO) 

'^h^.^h^^'u^Vd (DkT.j) - Da (ti^bc)) = -[Ka + La) T[^,^ + A.T^,,) - ^DdTb, - £(^bTc)a 

- i^^ad + (^ad) D'^Tbc — 2r^j '^£c)di^La^ + h^^bTc)' ^d 

--^'laibJ^cjd't^ +-^^{b-'-c)a ^^{bJ-c)a 

+2 {Qad + ^ad) K{bT,^ - 2 {^ad + ad) A^^T^^ 
-2 {Qa{b + 'Sa{b) Tc)dK''' + 2 {uJa,(b + O-a(b) 

2g) 2/5^ 2k^ 

+ ^/ia(6T,),A'' + — /l„(br,) '^gd - — g(6T,), . (All) 

These results apply for arbitrary (not necessarily small) scalar, vector and tensor fields. 



Appendix B: Infinitesimal frame transformations 

An infinitesimal frame transformation from the diad (u", n") to the diad (u", n") can be defined as a generalization 
of the corresponding general relativistic procedure [1^ as: 

Ua = Ua+Va + VUa , with M^Ua = n'^Va = , (Bl) 

na = na + la + mUa , with u°-la = n°-la = , (B2) 
where Ua, la, m are all O (1).'' The new diad also obeys 

vl^Ua = -1 , n'^ua = 1 , v^na = , (B3) 

which implies 

V ^ m . (B4) 

For Va — la = Q the above parameters define an infinitesimal orthogonal transformation (this is a 2-dimensional 
infinitesimal Lorentz boost). The parameters Va and la represent infinitesimal translations. These transformations 
represent gauge degrees of freedom, worth to explore in order to achieve particular tasks, for example to fulfill physical 
conditions, to conveniently close subsets of differential equations etc. 
The fundamental algebraic tensors hab and Sabc change accordingly: 

hab = hab + 2u^aVb) - '^n^alb) , (B5) 
Sabc = £abc - {ricEabe + ^n^aebjce) + {'^c^abd + 2u\^a'^b]cd) ■ (B6) 

The new 3-metric obeys hab'n"' — habu"' — 0. 



* The quantities denoted O (1) all vanish for identical transformation. Wc also assume a first order accuracy, thus all quantities O (2) = 
0(1)^ will be dropped. 
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The covariant derivatives of the new diad vectors can be invariantly decomposed as 

VaMb = -UaAb + KUaUb + KUaflb + UaKb + Lajlb + ^/iafa + ^ah + '^ab , (B7) 

Va^b = riaAb + KUaUb + KUaUb - UaKb + Lallb + ^hab + CJab + ^ab , (B8) 

which imply the following transformations for the kinematic quantities: 

K = K- Aav" + {Ka + La) t - TuK + m' , (B9) 

K = K- Aal" + (Ka - La^v" -mK -rh , (BIO) 

e = e + me + D'^L - ZMa - (^a - Ka) , (Bll) 

e = e + me + D''va+v^Aa-l''{Ka + La) , (B12) 

Ac = Ac + V^'AbUc + Klc - l^AbHc + + {oJac + <^ac) v" + m (Kc + Kc) + V(c) , (B13) 

Ac = Ac- l^AbUc - KVc + V^AbUc + + (Wac + ?ac) l" + m [Kc + Kc) + l[c) > (B14) 

Kc = Kc- Klc + UcV^Kb - Ucl^Kb + y^c + (o-ac + Wac) T + m {Ac + A^ + v\c) , (B15) 

Kc = Kc + KVc - net' Kb + UcV^Kb + jVc + {dac + Sac) + m (Ac + Ac) + /(c) , (B16) 

Lc = Lc- KVc - Klc + UcV'^La - Ucl^La - c 

+ ® Zc + Dctn - {deb + CJcb) + {aca + i^ca) l" , (B17) 

^cd = ^cd + mcjcd + D{cld) + [k^) - Ld)) - l(cAd) + 2v"-da(dUc) - "^l^^aidric) , (B18) 

^cd = (Ted + mdcd + D{cVd) - l{c (Kd) + Ld)) + vi^,-,Ad) - 2taa(dnc) + 2v"aa(d'"c) , (B19) 

Wed = oJcd + muicd + D[cld] + V[c (^Kd] - Ld]) - l[cAd] + 2v'^oja[dUc] - 2tuja[dnc] , (B20) 

Wed = Wed + mUJcd + L)[cVd] - l[c {Kd] + Ld]) + VlcAd] - 2rwa[d«-c] + 2z;"Wa[dWe] • (B21) 
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Similarly, the transformation laws for the gravito-electro-magnetic quantities are 

£ ^ £ + 2EJ'' + 2SJ'' , (B22) 

Ik ^ Sk- m£k - ^eik + Eav'^Uk - ej^nu + ekatn^v' + SkJ" ~ J'kaV^ , (B23) 

Ik ^ Sk- m£k + ^£vk - ial^'Uk + EaV-'Uk - EfcabH"/' + EkaV" - Fkal" , (B24) 

Uk ^ Uk- \ek ''''Eah - \ek "''£a^^^, - nkUJ" + UkUav" - "''SaVb + Ukav" - UkJ" , (B25) 

3 

= ^ki + 2u(k^i)ai}°' - 2n(^kJ^i)J"- - -£{kVi) 

+ lSi^kh) + miki + m£ki - Sahiknify' - Sabik'Hiff' , (B26) 

£ki = £ki + 2mTki + 2t;<fe£,) - Et^kk) + 2u(fc£,),w'' - 2n(fef ,),/'^ - 26,^^,.^;^ ''t , (B27) 

Ski = Ski + 2mTki - 2l(^k£i) + ^{kVi) + 2u(^kSi)av'' - 2n(^k£i)J'' - 2ea6(fc'H;)°w^ , (B28) 

nki = Uki + mUki + '^n^kh) - e(fc "''-^i)a4 + 2e(fe ''^£i)aVt 

'''SDaVb - 2n(fcHi)ar + 2«(feH,),u'^ , (B29) 
Uki = Uki + rnUki + \ni^kVi) + "'-F,),!;;, - 2e^J^ •'%)alb 

+£(, + 2u^kHi)aV'' - 2n(,H,)a/° . (B30) 

The matter variables transform as 

p = p - 2u'^g„ - 2mg , (B31) 

^ = + 2r^a - 2mq , (B32) 

P ^ P- \v''qa - \r^a , (B33) 

q = q + tqa-v''^a-m['p + l:) , (B34) 

5a = $a - PUa - g^a - 2pV a - 2tI abv" + UaV^qc - nj^qc - mTTa , (B35) 

TTa = Tfa - + 2pla ~ Tlla + 2TXabl^ - nJ^TTc + UaV'^TTc - niqa , (B36) 

TTab = T^ab + 2pV (a^b) " 2pl (^arib) " "^qiaVb) + '^T^ {a'^Ub)V d - 271 (^^"12^)^ - 2t: (Jb) ■ (B37) 



We have checked that in the particular case = = to, by applying the Lanczos equations (TfSll - lfT^ for eliminating 
the quantities O, La = —Ka and CTab, by suppressing the quantities related to the brane normal, in particular 
imposing Qab = 0, we recover the linearized form of the transformation laws for the kinematical and dynamical 
quantities (8, Oab, (jJab, ^a) and (p, p, qa, TTab) of Ref. [2^. Similarly, by employing Eqs. (|14p and ([T5|). we obtain 
the required transformations for the Weyl projections (Eab, Hab)- 

If we would like to apply the generic transformations derived in this Appendix in a brane- world scenario, we have to 
impose = = m on the brane (in order to preserve the vector field n° at y = 0, which defines the brane), however 
the derivatives along the off-brane direction (the derivatives denoted by prime) of these quantities can be different 
from zero even at ?; = 0. 

Appendix C: Gravitational evolution and constraint equations on an asymmetrically embedded brane 

The brane equations describing the gravitational dynamics are 



= e - D'^ka + + f j - "^K^Aa + daho''^ - , 

~2A,kad'''' + '^dabS"' + + l) "-^^^^ + ^-^f + ^-Kak'^ ' aabk'^k" 

x4^1 4£ ^ ^ ^ ^ ^ 

£{k) + 3®^*= ~ 3^'^^ S"^*^ ~ D^^ka - £ka^"' " (Wfca - ka) + ^ {ka)K°' 

-{k + ^ D'dkb + ^Vefca/ + '^(^k + ^^ Dke + ^(^k + 1 j kk - 2^»A^a + ^"£>a^fe 

+ i^fc£)»^„ - 2k''A^kka) - UbadkK" + ^kkUabd"-" + a^'iPfeaab - Ub'^D^'dka - \dk"DaQ 

+ecabK''J'dk + efe "^ifcWa^f, " - ^^r'^fe) + y^fc (p + Stt - 3p) + — I if - - 1 TTfe 

- —9 \2Kk + /Tfe j + — TTfeaA" - — (tT - p) Afc - — CTfcaTr" , 

B2 - - 1~ 

e - + _ + eiv: - AM„ - 2w„a;'' + uabo"'' - K'^Ka - £ - -K+ — {p + n +p) , 

1 20 



29 i/^e\^ 1^^ 

1 1 _ 
+W(oWb) + Uc^a'^i,)" + -j^cia^b)'' + -^afc + 2'^"'' S"''""^ ' 

= D^Wa - A^W" , 

= -D(cWfc) + eab{kD^(Jc)°' + 2A(cWfc) + Hab , 

2 20 ^ ^ 
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= -B(fej) - -S{kj) - SabikD^Hj^ + -£'(fc^j) + QEkj - —£kj + ^ik^j) - -^^''^ 

1^^3^^^ '^1 ^ 

--^CrjkK°'Ka + -Kf^k^j-jaK"- + K(k'^j)aA°' + da{jAk)K'' - i^i^a j)a — ^T^ikj) 

~2 ~2 ^ ~2 ~2 ~2 
3--D(fcgj> + y7f(fe-ftrj) - -^g(/c^j) - y (P + P)<yjk - -g-©7rjfc - y^^o" ('^fc>a + (^fc>a) , (dl) 

= + ^ah(kD°'E-'^ + -eab{kD"£j^ + QHkj - 3aa{kHj^ - uJa{kHj) 

^ - ) ^abikD'^dj^ + -eab{kD''dj)cd''^ + -eab{k'<fj}cD"'a'''' 



+ea;(fei?,) + l^^, ""'a.^aKb - ^e^bikK.^D^K'' - ^e^.^.K^'D'^K^^ 

+ \^{k °'^'^3)b'^a''Kc - ^d(°'UJk)Ka - y £a6(fe-D°7r5) - K^ji^k) " ye(fe °'^(^j)aQb , (C12) 



1^1 B 1 2B 

= D-E,k - ^D-£,k + ^Dk£ - SHkaCo'' + '^'^H^.at + j^k - ^ {ScJka + cTka) - -^Dk® 

-IdakD'^ - l) - ^ - l) - l^'D'^^ck + l^''''Dkaab - l^kbD'^dJ' 

-^K-DkKa + iKaD-Kk + ^KkD'^Ka + ^Sk ""KeLo,^/ - ^dj'k'^akb - \^QKk + ^<JkbK' 

+@e^,''dK,Wd + -dk''Ka + ^D'^'Hak - -^DkCp-li +P) + —Qqk - {'i^ka + <Tka) , (C13) 

= D'^Hak + \ek """DaSb - y ^fe + ■iEkaio'^ - £fe "'-E„,a, + if afcw" + \ek """Sac^b " - \K-'^cKk 

'^'KbD^e - l^fe '^''DSb - IsabkK-D'K'' - '^^'a^'DbK, + ^ + |^ u;, - ]^k^k-oJk 

+ \h "'^dadc "^b ' - y ^feaW" + y Efc "'D.qb + ^ (p + p) Wfe - '"'^a^Obc • (CM) 

We note that this system of equations is vahd on both sides of the brane. 
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Appendix D: Kinematical, gravito-electro-magnetic and matter variables for Bianchi I brane-worlds 

In this Appendix we rewrite the Bianchi I brane-world, presented in Ref. [49] in terms of our variables. The 
brane- world solution contains an unspecified function V (y) of a Gauss normal coordinate y (which is however not 
related to the brane normal). 

The kinematical quantities appearing on and outside the brane are presented in Tables |TT] and IIIIl respectively, 
while the gravito-electro-magnetic quantities on and outside the brane are given in Tables IIVI and [V] Tables IIIIl and 
|V] contain quantities, which were not computed in Ref. [A^ . 

TABLE II: Brane kinematical quantities (first column) for the brane-world ^] (second column; notations and Eq. numbers 
are from this reference). 



kinematic quantity 


for the metric in Ref. [49] 







e 


Q, as given by Eq. (63) 


O'ab 


OAB, as given by Eqs. (64)-(65) 








TABLE III: Off-brane kinematical quantities (first column) for the brane-world [1^ (second column; notations are from this 
reference) . 



kinematic quantity 


for the metric in Ref. [49] 


K 


e / 1 Co 1 VV' \ 


K 


1 / Vu' 1 VCo 1 V' \ 


e 


e / 3u' Co 


Aa 





Ka 





Ka 





La 





2a 





'Sab 


3 

Y^aiCiaCit , (Ti = / „ {Co + 3Ci) 



TABLE IV: Brane gravito-electro-magnetic quantities (first column) for the brane-world [i^ (second column; notations and 
Eq. numbers are from this reference). 



Weyl quantity 


for the metric in Ref. [49] 


£ 


-K^W given by Eq. (53) 


£a 





£ab 


-k^Vab given by Eqs. (54)-(55) 



The notations for the brane matter variables (after the straightforward change in the indices from lower case to 
upper case letters) are identical in this paper and in Ref. [!§], with qa = 0. There are no off-brane matter variables. 
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TABLE V: OfF-brane gravito-electro-magnetic quantities (first column) for tfie brane-world [i^ (second column; notations are 
from this reference). 



Weyl quantity 


for the metric in Ref. [49| 


£a 










£ab 


E E,erae,k , E, = [(Co + SC.) u' ~ - 2 (Cg + 3Cf)] 

i = l 

+ ,^.(ly2^ [iCo + SCO u'-C + 4^ (Co - CO] 


J'ab 




'H.ab 





Hab 
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